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An  analysis  of  water  management  in  air  breathing  small  polymer  electrolyte  membrane  fuel  cells  (PEM- 
FCs)  is  presented.  Comprehensive  understanding  of  flooding  and  dry-out  limiting  phenomena  is  presented 
through  a  combination  of  analytical  modeling  and  experimental  investigations  using  a  small  PEMFC  pro¬ 
totype.  Configurations  of  the  fuel  cell  with  different  heat  and  mass  transfer  properties  are  experimentally 
evaluated  to  assess  the  impact  of  thermal  resistance  and  mass  transport  resistance  on  water  balance. 
Manifestation  of  dry-out  and  flooding  problems,  as  limiting  phenomena,  are  explained  through  a  ratio 
between  these  two  resistances.  Main  conclusions  are  that  decreasing  the  ratio  between  thermal  and  mass 
transport  resistance  under  a  certain  point  leads  to  flooding  problems  in  air  breathing  PEMFC.  Increasing 
this  ratio  leads  to  dry-out  of  the  polymer  electrolyte  membrane.  However,  too  high  thermal  resistance 
or  too  low  mass  transport  resistance  reduces  the  limiting  current  by  pushing  forward  the  dry-out  prob¬ 
lem.  This  work  provides  a  framework  to  achieve  the  proper  balance  between  thermal  rejection  and  mass 
transport  to  optimize  the  maximum  current  density  of  free  convection  fuel  cells. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

One  of  the  main  problems  in  the  development  and  commercial¬ 
isation  of  micro  polymer  electrolyte  membrane  fuel  cells  (PEMFCs) 
as  small  power  sources  (<10W)  resides  in  the  miniaturisation  of 
auxiliary  components  that  must  be  scaled  down  with  the  stack. 
External  subsystems,  such  as  the  blower,  humidification  unit  and 
heat  exchanger,  increase  the  complexity  for  balance  of  the  plant  of 
the  fuel  cell  system  and  typically  have  low  efficiencies  and  power 
density  when  miniaturised.  This  suggests  that  micro-PEMFC  [1-4] 
may  not  lead  to  viable  power  sources  unless  research  is  done  on 
miniature  balance  of  plant  components  or  passive  approaches  are 
developed. 

Running  PEMFCs  in  free  convection  mode  significantly  reduces 
the  number  of  subsystems  required.  Absence  of  direct  control  over 
the  air  stream  in  the  PEMFC  makes  temperature  and  humidity 
control  an  important  problem.  In  fact,  air  breathing  fuel  cells  consis¬ 
tently  achieve  lower  performance  than  cells  operated  under  forced 
convection  mode.  Many  studies  have  been  conducted  on  fuel  cells 
running  in  free  convection  mode,  either  through  modeling  or  exper¬ 
imentation  [5-19].  Main  conclusions  are  that  in  this  operation 
mode,  problems  limiting  performances  are  due  to  mass  transport 
and  water  management.  Also,  the  importance  of  the  surrounding 
conditions  is  highlighted  [7,8,13]. 
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Mennola  et  al.  [14]  concluded  that  for  cell  temperatures  below 
40  °C,  the  limiting  problem  was  the  system’s  inability  to  remove 
water.  Modroukas  et  al.  [15,16]  conducted  a  study  of  air  breath¬ 
ing  PEMFC  running  at  low  temperature  (20-30  °C).  Under  these 
operating  conditions,  the  limiting  factor  was  also  water  accu¬ 
mulation  which  caused  flooding  of  the  cell.  To  increase  cell 
performance,  utilisation  of  a  passive  water  control  system  based 
on  hydrophilic  material  that  enhances  water  removal  from  the 
cell  was  demonstrated.  O’Hayre  et  al.  [17]  ran  a  free  convection 
cell  under  temperatures  ranging  from  10  to  80  °C  and  demon¬ 
strated  the  dependence  of  performance  on  humidification  of  the 
cell.  At  temperatures  higher  than  70  °C,  membrane  dry-out  was 
found  to  be  the  limiting  phenomenon.  Morner  and  Klein  [18]  also 
concluded  that  membrane  humidity  is  one  of  the  most  significant 
factors  in  determining  fuel  cell  performances  in  free  convection 
mode.  Jeong  et  al.  [8]  also  concluded  that  the  performances  of 
air  breathing  fuel  cells  are  strongly  dependant  on  water  manage¬ 
ment. 

This  past  work  highlighted  two  interesting  points.  First,  water 
management  problems  seem  to  be  the  limiting  phenomena,  either 
by  flooding  the  cathode  domain  or  by  dry-out  of  the  electrolyte.  Sec¬ 
ond,  the  operating  temperature  is  important  in  defining  the  limiting 
problem.  Running  in  free  convection  mode  implies  a  strong  depen¬ 
dency  between  heat  transfer,  mass  transfer  and  temperature  of  the 
system.  Packaging  dimensions  and  material  properties  are  then  of 
considerable  importance,  especially  for  fuel  cell  systems  without 
air  flow  and  temperature  control,  such  as  portable  air  breathing 
products. 
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Nomenclature 

Cp  specific  heat  (J  kg-1  K-1 ) 

D  diffusion  coefficient  (m2  s-1 ) 

E0cv  reference  reversible  fuel  cell  voltage 

F  Faraday’s  constant  (96,487  C  mol-1 ) 

g  gravity  constant  (9.8  m  s-2 ) 

Agf  Gibbs  free  energy  change 

h  free  convection  heat  or  mass  transfer  coefficient 

(W  m-2  K-1 )  or  (m  s-1 ) 

Hfg  heat  of  water  condensation  (40  kj  mol-1 ) 

I  current  density  (A  m-2 ) 

/0  exchange  current  density  (A  nrr 2 ) 

I\  limiting  current  density  (A  m-2 ) 

/i_ max  maximum  limiting  current  density  (A  m-2 ) 

/i_ ideal  ideal  limiting  current  density  (A  nrr 2 ) 
k  thermal  conductivity  (W  nr1 I<-1 ) 

L  characteristic  length  for  heat  and  mass  transfer 

fn  mass  flow  (kgs-1) 

M  molar  weight  (g  mol-1) 

N  molar  flux  (mol  s-1  m-2 ) 

Q  heat  (W) 

Qg  heat  generated  by  fuel  cell  (W) 

R  universal  gas  constant  (8.314J  mol-1 1<) 

^electric  fuel  cell  electrical  resistance  (£2) 
iconic  fuel  cell  ionic  resistance  (£2) 

RH  relative  humidity 

S  surface  (m2) 

T  temperature  (K) 

w  mass  fraction 

Wqdl  GDL  thickness  (m) 

x  molar  fraction 

Greek  letters 

ace\\  charge  transfer  coefficient 

^membrane  Nation®  membrane  thickness  (m) 

£  porosity 

cr rh  membrane  conductivity  (£2-1  m-1 ) 

Subscripts 

a  active  area 

c  cathode 

cell  fuel  cell 

cr  critical 

GDL  gas  diffusion  layer 

m  mass  transport 

packaging  cell  subsystems  excluding  cathode  and  GDL 
s  GDL  interface  with  environment 

T  thermal 

a  species:  N2,  02,  H20 

(3  species:  N2,  02,  H20 

oo  environment 


In  this  paper,  we  develop  a  simple  one-dimensional  model  of  a 
free  convection  PEMFC  with  a  planar  configuration  that  integrates 
heat  and  mass  transfer.  This  model  will  be  validated  experimentally 
and  then  used  along  with  other  experimental  results  to  analyse 
the  impact  of  thermal  resistance  and  mass  transport  resistance.  To 
do  so  various  configurations  of  the  prototype  with  different  heat 
and  mass  transfer  properties  will  be  evaluated  experimentally.  The 
manifestation  of  dry-out  and  flooding  as  limiting  phenomena  will 
be  presented  in  term  of  ratio  between  heat  and  mass  transfer.  This 
work  unifies  diverging  conclusions  from  previous  work  [14-17]  on 
the  manifestation  of  the  limiting  problems.  Also,  the  importance 


of  the  system’s  packaging  will  be  highlighted  and  experimentally 
shown. 


2.  Free  convection  planar  PEMFC  model 

2.1.  Free  convection  heat  and  mass  transfer  coefficient 


Free  convection  on  exposed  surfaces  of  a  planar  PEMFC  is 
characterized  by  heat  and  mass  transfer  resulting  from  the  elec¬ 
trochemical  reaction  in  the  cell.  Oxygen  consumption  and  water 
production,  resulting  from  cell  operation,  will  generate  concentra¬ 
tions  gradients  and  density  changes  in  the  air.  Fleat  loss  from  the 
cell  creates  a  temperature  gradient  which  generates  a  density  varia¬ 
tion  in  the  air  surrounding  the  cell  and  creates  convection.  All  these 
gradients  will  cause  heat  and  mass  transfer  between  the  cell  and 
its  environment.  An  exact  free  convection  model  should  consider 
both  heat  and  mass  transfer  phenomena  combined.  But,  as  a  first 
approximation,  heat  transfer  and  mass  transfer  will  be  considered 
separately,  related  by  the  heat  and  mass  transfer  analogy.  In  this 
section,  the  heat  transfer  coefficient  (hj)  and  air  species  mass  trans¬ 
fer  coefficient  (h“)  will  therefore  be  determined  using  empirical 
correlations. 

Natural  convection  heat  transfer  can  be  described  by  the  non- 
dimensional  heat  transfer  coefficient,  referred  to  as  the  Nusselt 
number  ( NuL ).  For  natural  convection  driven  by  heat  transfer  only, 
empirical  correlations  depending  on  geometries  and  boundary  con¬ 
ditions  express  the  Nusselt  number  as  a  function  of  the  Grashof 
(GrL)  and  Prandtl  (Pr)  non-dimensional  numbers  [20]: 


Nul  =  - I^"=f(CrL,Pr)  (1) 

Free  convection  from  a  planar  PEMFC  will  simply  be  represented 
as  a  vertical  heated  plate.  For  this  situation,  the  heat  transfer  cor¬ 
relation  is  [20]: 

Nul  =  0.54( CrL  Pr)1/4,  for  104  <  CrL  Pr  <  107  (2) 


GrL  Pr  =  RaL 


mjs  -  Too)!3 

va 


(3) 


As  presented  in  Eqs.  (2)  and  (3),  only  the  temperature  gradient  is 
considered  to  lead  to  a  density  change  in  the  air.  Species  concen¬ 
tration  variation  will  also  induce  density  gradient,  but  this  effect  is 
not  considered  here  for  simplicity.  As  presented  by  O’FIayre  et  al. 
[17],  the  thermal  gradient  is  the  main  factor  in  natural  convection 
of  air  breathing  PEMFC  and  accounts  for  about  75%  of  the  density 
gradient.  This  approximation  leads  to  an  underestimated  value  of 
the  heat  transfer  coefficient. 

The  Sherwood  (Sh^)  non-dimensional  number  can  be  used 
to  describe  free  convection  mass  transfer.  This  number  can  be 
empirically  correlated,  depending  on  geometries  and  boundary 
conditions,  with  the  non-dimensional  Grashof  (Gri.)  and  Schmidt 
(Sc)  numbers: 


Shi  =  feirzE  =/(CrL,  Sca)  (4) 

In  order  to  use  the  analogy  between  heat  and  mass  transfer  for 
a  fluid  flows  where  there  is  simultaneous  heat  and  mass  transfer, 
the  Lewis  (Le)  number,  which  is  the  ratio  between  heat  and  mass 
diffusivity,  should  be  close  to  unity.  For  diffusion  of  oxygen  and 
water  vapour  in  ambient  air  at  25  °C,  the  Lewis  number  is 

Leo2  =  ^  =  l-04  (5) 

LeH2o  =  =  0.82  (6) 
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Considering  that  the  Lewis  number  values  are  near  unity,  the 
analogy  between  heat  and  mass  transfer  will  be  used  here.  An 
equation  for  each  air  species  (N2,  02,  H20)  can  then  be  obtained 

Sh?  =  0.54(Gr“SCc,)1/4,  for  104  <  Cr“  Sca  <  107  (7) 


The  term  GrL  Sca  from  Eq.  (7)  is  calculated  using  a  relation  pro¬ 
posed  by  Li  et  al.  [9]: 


Gr*Sca  =  Ra<l  = 
where 

M3ir 
Ya 


gKa(Wa-S  -  Wq-cJI3 

—n 


waMair  +  ((l  -  wa)Ma ) 


(8) 


(9) 


Knowing  the  mean  value  of  Nul  and  ShL,  it  is  possible  to  calculate 
the  mean  value  of  the  heat  and  mass  transfer  coefficients: 


hj  - 


h«  = 


0.5  4(GrLPr)V4kai[ 

L 

0.54  (CrfSCc,)1^^ 


(10) 

(11) 


2.2.  Model  scope  and  assumptions 


The  scope  of  the  model  is  to  identify  the  manifestation  of  lim¬ 
iting  problems  at  the  cathode  of  an  air  breathing  PEMFC;  either 
water  management  problems  or  mass  transport  limitations.  Quali¬ 
tative  results  that  show  dominant  physics  are  expected,  as  opposed 
to  precise  predictive  results.  This  model  is  inspired  of  those  pre¬ 
sented  by  Refs.  [11,12,14,15,17,19].  The  main  assumptions  used  to 
develop  the  model  are 


for  membrane  water  content  and  its  related  ionic  conductivity. 
Water  absorption,  desorption  and  transport  in  the  Nation®  mem¬ 
brane  [21,22]  are  not  consider  here.  Future  work  should  consider 
water  transport  in  the  entire  membrane  electrode  assembly  (MEA), 
including  water  balance  at  the  anode. 

2.3.  Model  domains 

A  control  volume  approach  is  used  to  consider  energy  and  mass 
balance.  The  model  is  divided  in  three  domains.  The  first  domain  is 
the  GDL  which  goes  from  cathode’s  catalyst  region  to  the  interface 
of  the  GDL  with  the  environment.  At  the  interface  of  the  catalyst 
with  the  GDL,  the  boundary  conditions  represent  the  electrochem¬ 
ical  reactions  occurring  in  the  cell,  namely:  heat  generation  (Qg), 
water  production  (Nh2o)  and  oxygen  consumption  (Nq2).  These 
fluxes  are  directly  related  to  the  cell  current  density  (/).  The  sec¬ 
ond  domain  is  the  GDL  interface  with  the  environment.  Boundary 
conditions  for  this  domain  are  environment  species  molar  fraction 
(*oo-a)  and  temperature  (Too).  The  third  domain  represents  rest  of 
the  fuel  cell,  such  as  the  packaging,  that  can  exchange  heat  with 
the  environment  at  Too;  mass  transport  only  occurs  in  domains  1 
and  2.  The  presence  of  this  third  domain  is  important  considering 
the  analysis  presented  here,  which  includes  the  packaging  of  the 
system  and  not  only  the  cathode  and  the  GDL.  The  diagram  in  Fig.  1 
presents  heat  and  mass  transfer  fluxes  and  sketches  the  trends  of 
species  molar  fraction  and  temperature  within  the  domains. 

In  the  GDL,  heat  transfer  is  modeled  as  conduction  only  (Qc). 
Mass  transfer  occurs  via  species  diffusion  through  the  porous 
media.  The  governing  equation  for  the  displacement  of  oxygen, 
water  vapour  and  nitrogen  is  given  by  the  Maxwell-Stefan  diffusion 
equation  for  multi-component  systems  [23]. 


•  one-dimensional; 

•  steady  state; 

•  dry  hydrogen  at  the  anode; 

•  no  water  crossing  the  membrane  to  the  anode; 

•  uniform  activity  over  the  cathode  surface; 

•  infinitely  thin  cathode; 

•  driving  force  for  mass  transport  in  the  gas  diffusion  layer  is  by 
diffusion  only  [15]; 

•  uniform  water  content  across  the  Nation®  membrane; 

•  equilibrium  between  relative  humidity  at  the  cathode  and  water 
activity  in  the  membrane. 

An  important  assumption  of  the  model  is  the  absence  of  satura¬ 
tion  effects  due  to  water  condensation  within  the  gas  diffusion  layer 
and  in  the  catalyst  region.  Modroukas  [15]  showed  the  importance 
of  this  aspect  to  model  a  free  convection  PEMFC  running  under  free 
convection  at  low  temperature.  Under  these  conditions  water  accu¬ 
mulation  is  the  dominant  problem  that  causes  cell  failure.  Since 
the  focus  here  is  to  identify  if  flooding  will  occur  and  not  to  predict 
the  performance  once  it  has,  the  model  does  not  include  the  effect 
of  water  accumulation.  Therefore,  the  model  will  not  include  the 
impact  of  relative  humidity  (RH)  going  over  100%  and  subsequent 
water  formation  that  would  impede  oxygen  diffusion  to  the  cata¬ 
lyst  region.  The  important  aspect  that  will  be  captured  is  when  the 
water  partial  pressure  is  higher  than  the  water  saturation  pressure, 
to  identify  when  water  will  condense  in  the  cathode  and  flooding 
will  occur.  To  appropriately  do  so,  heat  of  water  condensation  (Hfg) 
will  be  taken  in  account  as  heat  generated  by  the  cell. 

Another  impact  of  the  assumptions  presented  above  is  related  to 
water  content  in  the  Nation®  membrane.  The  last  two  assumptions 
listed  above  imply  that  water  conditions  within  the  membrane  will 
be  fixed  by  the  temperature  and  partial  pressure  of  water  at  the 
interface  between  the  cathode  and  the  GDL.  This  assumption  limits 
the  predictive  accuracy  of  the  model  but  captures  the  tendency 


Qc  =  -k  gdlSc 


E 


Ts-Tc 

Wgdl 
RTC 

™(a-P)_eff 


(xpNa  —  xaNp)  —  —  Vxa, 


P=o2,h2o,n2 

fora  =  02,  H20,  N2 


(12) 


(13) 


Binary  diffusion  coefficients  of  species  (Da_p)  are  calculated 
using  relations  from  Bird  [23].  In  order  to  take  into  account  porous 
media  effects,  the  Bruggemann  correction  factor  is  used  to  correct 
the  binary  diffusion  coefficients  for  GDL  porosity  (£Gdl)  [14]. 

D“-p = ap  (70^)  (PcraPcrp)1/3(:rcraTcrp)5/,2\/i+i 

(14) 

3 /2 

^(a-(3)_eff  =  ^GDL^a-P 

On  the  exposed  surface  of  the  GDL,  heat  and  mass  transfer  are 
driven  by  free  convection.  Mean  values  of  heat  and  mass  transfer 
coefficients  (hj  and  )  are  used.  The  following  are  the  governing 
equations  that  link  heat  rejection  and  species  movement  on  the 
GDL  surface: 

Qc  —  hjSa{Ts  —  T^)  (16) 

Y  (xpN«  -  x„Np)  =  -AxaT^.fi“  (17) 

P=o2,h2o,n2 

An  important  value  pictured  on  Fig.  1  is  QoUt,  which  accounts  for 
all  other  heat  loss  mechanisms  in  the  system  except  that  from  the 
GDL  surface  (Qc).  A  significant  part  («90%)  of  the  heat  generated 
by  the  fuel  cells  will  be  rejected  in  the  environment  through  the 
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(Free  Convection)  (Free  Convection) 


Fig.  1.  Domains  for  the  free  convection  planar  PEMFC  model  with  heat  transfer  (top)  and  mass  transfer  (bottom)  phenomena.  The  free  convection  transfer  processes  at  the 
cathode  are  depicted  on  the  right  hand  side,  while  the  other  transfer  processes  for  the  anode  and  the  rest  of  the  packaging  are  shown  on  the  left. 


rest  of  the  system  (packaging)  in  contact  with  the  environment  and 
through  heating  of  the  system’s  hydrogen  supply  (Qout). 

Qout  —  C^cell  —  5a)/li_packaging(  ^packaging  —  ^oo) 

+  mH2Cp(Tc-T00)  (18) 


and  the  material’s  internal  resistance.  Ionic  resistance  is  mainly 
related  to  the  electrolyte  membrane  thickness,  water  content  and 
temperature.  The  ionic  conductivity  dependency  with  water  con¬ 
tent  of  the  membrane  will  be  estimate  using  an  empirical  relation 
developed  by  Sone  et  al.  [25],  which  relates  ionic  conductivity  of 
Nation®  with  relative  humidity. 


2.4.  Relation  between  heat  and  mass  transfer  and  the 
electrochemical  reactions 

A  fundamental  part  of  the  model  is  to  gather  heat  and  mass 
transfer  physics  with  the  electrochemical  reaction  that  occurs 
within  the  cell.  Environment  conditions  combined  with  heat 
generated  and  mass  fluxes  are  the  input  necessary  to  establish 
temperature  and  species  partial  pressure  of  the  cell  for  an  oper¬ 
ation  point.  These  last  two  inputs  are  driven  by  the  electrochemical 
reaction  of  the  cell  [24]. 

Using  the  stoichiometry  of  the  reaction,  it  is  possible  to  establish 
the  following  molar  flux  for  oxygen  consumption  and  water  vapour 


generation. 

(19) 

M  1 

Nh2o  =  2P 

(20) 

In  order  to  describe  the  fuel  cell  polarization  curve  (V-J  curve), 
voltage  drop  must  be  determined.  There  are  three  main  sources  of 
voltage  drop:  reaction  kinetics  (%),  Ohm’s  law  (rjlR)  and  mass  trans¬ 
port  (?7mt)-  Voltage  drop  from  reaction  kinetics  can  be  described  by 
the  Tafel  equation  [24]: 

where  J0  is  the  exchange  current  density  and  acell  is  the  charge 
transfer  coefficient.  These  values  will  be  fitting  parameters  in  the 
model.  Voltage  drop  from  Ohm’s  law  is  due  to  internal  electrical 
resistance  and  ionic  resistance  of  the  Nation®  membrane.  Electrical 
resistance  is  mainly  due  to  contact  resistance  between  components 


4lR  —  (^electric  +  ^ionique) 

p  _  ^membrane 

SatfRH 

ctrh  =  7.46x3  -  7.45x2  +  3.13x  -  0.378 


(22) 

(23) 

(24) 


where  x  is  the  relative  humidity.  As  presented  in  the  assumptions, 
this  value  will  be  approximate  by  relative  humidity  at  the  cath¬ 
ode  (RHC)  because  equilibrium  is  considered  between  the  relative 
humidity  at  the  cathode  and  the  water  activity  in  the  membrane. 
Also,  Eq.  (24)  is  accurate  only  for  a  temperature  of  30  °C.  As  pre¬ 
sented  in  Ref.  [25]  ionic  conductivity  of  a  Nation®  membrane  not 
only  depends  on  its  configuration  and  relative  humidity  but  also  on 
its  temperature. 

RHc  =  ^c~H2°-(gas)P  (25) 

‘  saturation,  H20 

Using  saturation  pressure  values  from  Van  Wylen  et  al.  [26],  a 
third  order  polynomial  function  was  defined  to  fit  saturation  pres¬ 
sures  between  15  and  70  °C. 

^saturation, h20  =  0.164TC3  +  145. 327)?  +  43023. 33TC 

-4263759.49  (26) 


Voltage  drop  due  to  mass  transport  comes  from  the  decrease  in 
oxygen  partial  pressure  at  the  catalyst  surface  [24]. 


4mt  = 


-RTC 
4  F 


In 


f  *o2-c 

y  Xq2-oo 


(27) 
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It  is  possible  to  write  the  polarization  value  of  the  fuel  cell  as 
follows  [11]: 


V  =  Eocv  -  Pk  -  P\r  -  Ptm 

where  Eocv  is  the  fuel  cell  open  circuit  voltage  [24]. 
-Agf 


Eo  = 


2  F 


(28) 

(29) 


System  loss  will  result  in  heat  generation  (Qg)  proportional  to 
voltage  drop  (rjK,  r\ iR  and  tjtm)-  We  consider  that  water  is  pro¬ 
duced  in  vapour  form,  so  to  be  thermally  consistent;  we  calculate 
Eq  using  lower  heating  value  (LHV)  of  the  Gibbs  free  energy  varia¬ 
tion  (Agf).  Heat  released  during  water  condensation  inside  the  cell 
is  accounted  for  using  heat  of  water  condensation  (Hfg).  The  heat 
generation  rate  can  be  expressed  as 


Qc  =  /Sa 


^3fg*c-H20_liquid  \ 

2 F  J 


(30) 


If  RHC  <100%,  then  the  molar  fraction  of  liquid  water  at  the  cath¬ 
ode  is  xc_H2o_iiquid  =  0,  otherwise  xc-H2o_iiquid  is  calculated  using 
saturated  value  of  water  molar  fraction  (xc_h2o  (saturated))1 

*c-H20_liquid  =  *c-H20  —  *c-H20  (saturated)  (31) 

Using  the  equations  developed  above,  an  iterative  approach  can 
be  used  to  find  a  temperature  value  of  the  cell  (Tc)  that  will  match 
heat  generated  through  electrochemical  reactions  (Qg)  with  heat 
loss  to  the  environment  (Qc  and  QoUt).  To  do  so,  an  approximate 
value  of  Tc  will  be  fixed  and  used  to  calculate  the  temperature  at  the 
surface  of  the  GDL  (Ts)  based  on  conduction  through  the  GDL  and 
the  stainless  steel  grid  that  are  over  the  cathode.  The  temperature  of 
the  rest  of  the  packaging  (Tpackaging)  is  calculated  as  heat  conduc¬ 
tion  through  a  typical  thickness  of  the  fuel  cell  packaging.  These 
temperatures  will  then  be  used  to  calculate  heat  transfer  coeffi¬ 
cients  hj  and  hT  packaging  with  Eq.  (10).  Qc  and  QoUt  will  then  be 
determined  using  Eqs.  (16)  and  (18).  It  should  be  noted  that  this 
simple  modeling  of  heat  loss  from  the  system  is  corrected  by  a  fit¬ 
ting  factor  (Cheat)  that  is  determined  from  the  experimental  results 
of  the  baseline  prototype,  and  used  throughout.  The  energy  balance 
presented  in  Eq.  (32)  is  done  at  a  specific  operating  point,  defined 
by  the  current  density  (/)  which  also  fixes  oxygen  and  water  vapour 
partial  pressure  by  imposing  the  corresponding  molar  fluxes  as  pre¬ 
sented  in  Eqs.  (19)  and  (20).  Knowing  species  partial  pressure  and 
an  approximation  of  temperatures  within  the  system,  voltage  loss 
(V)  and  heat  generated  (Qg)  are  calculated  using  Eqs.  (28)  and  (30). 
Iterations  will  be  done  around  Tc  as  long  as  the  energy  balance,  Eq. 
(32),  is  not  satisfied. 

Qg  —  Qc  +  Qout  (32) 


3.  Experimental 

3.1.  Small  PEMFC  prototype 

In  order  to  validate  the  model,  a  small  PEMFC  was  designed  and 
fabricated.  Fig.  2  presents  an  exploded  view  of  the  prototype  with  a 
description  of  the  various  components.  External  dimensions  of  the 
prototype  are  3.5  cm  long  by  1.9  cm  large  and  0.635  cm  thick.  Active 
area  of  the  cell  is  2  cm2.  Table  1  summarises  the  values  used  in  the 
analytical  simulation  of  the  prototype. 

3.2.  Experimental  set-up  and  results 

Tests  were  run  using  a  fuel  cell  test  station  fabricated  by  Fuel  Cell 
Technologies,  Inc.  Dry  hydrogen  was  supplied  to  the  anode  using  a 
stoichiometry  ratio  of  3  and  atmospheric  air  was  supplied  to  the 
cathode  by  free  convection  with  the  ambient.  During  testing,  the 


(a) 

(b) 

(c) 


Fig.  2.  Exploded  view  of  the  small  PEMFC  prototype,  (a)  Stainless  steel  support  plate 
for  cathode,  (b)  Grafoil®  sheet,  (c)  stainless  steel  grid  for  cathode,  (d)  GDL  (SGL 
Carbon  10  BB),  (e)  MEA:  catalyst  coated  NAFION®  NRE  212  and  loading  of  0.5  mg 
Ptcm-2  from  Ion  Power,  Inc.,  (f)  seal  to  prevent  lateral  diffusion  of  hydrogen  in  anode 
GDL,  (g)  stainless  steel  anode  channel,  (h)  seal  for  anode  flow,  and  (i)  Plexiglass 
support  plate  for  anode. 


polarization  curve  ( V-I )  was  measured  as  well  as  the  temperature 
and  relative  humidity  of  the  cathode  GDL’s  exposed  surface.  It  is 
important  that  these  measurements  be  made  as  close  as  possible  to 
the  GDL  but  geometrical  limitations  did  not  allow  measurements 
to  be  taken  directly  on  GDL’s  top  surface.  Hence,  we  expect  the 
actual  GDL  temperature  to  be  slightly  higher  than  the  measure¬ 
ment,  leading  to  some  discrepancies  between  results  and  modeling. 
Measurement  uncertainty  in  current  density  and  voltage  are  small, 
respectively  ±0.003  A  and  ±0.003  V.  For  temperature  and  relative 
humidity,  we  evaluate  the  uncertainty  to  be  respectively  less  than 
±1  °C  and  ±7.5%  RH,  considering  sensor  accuracy  and  measurement 
repeatability. 

Experimental  results  were  taken  in  steady-state  conditions 
except  for  points  having  a  relative  humidity  of  100%.  These  operat¬ 
ing  points  presented  water  accumulation,  which  mean  that  water 
flux  generated  by  the  reaction  was  higher  than  the  maximum  water 


Table  1 

Model  input  values 


Model  parameter 

RHoo 

As  specified 

Too 

As  specified 

^electric 

0.2  £2 

Eocv 

1.185  V 

wGDl 

0.38  mm 

£gdl 

0.84 

Sa 

2  cm2 

Scell 

20  cm2 

L 

3.5  cm 

^membrane 

70|jim 

<^cell 

0.22  [24] 

/<GDL 

lOWm-1  K-1  [17] 

Fitting  parameter 

Io 

0.07  Am-2 

Cheat 

2.17 
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Fig.  3.  Comparison  between  experimental  and  simulation  results  for  free  convection  fuel  cell  operated  at  ambient  air  of  23.5  °C  and  40%  RH  as  a  function  of  current  density: 
(a)  polarization  curves,  (b)  relative  humidity  on  surface  of  cathode  GDL,  (c)  fuel  cell  membrane  resistance  calculated  by  Eqs.  (23)  and  (24),  based  on  relative  humidity  at  the 
surface  of  the  cathode  GDL  for  the  experimental  results,  and  (d)  temperature  on  surface  of  cathode  GDL 


flux  removed  by  free  convection.  Therefore,  long-term  operation 
under  these  conditions  would  create  water  accumulation  at  the 
cathode,  significantly  reducing  oxygen  diffusion  and  limiting  the 
achievable  current  density  due  to  mass  transport  resistance.  Results 
under  these  conditions  were  taken  before  water  accumulation,  so 
these  results  are  not  achievable  in  long  term,  steady-state  opera¬ 
tion. 

Fig.  3  presents  experimental  results  for  the  small  PEMFC  pro¬ 
totype  compared  to  results  obtained  with  the  model  developed 
previously.  The  prototype  limiting  current  density  (I\)  occurs  at 
360  mA  cm-2  and  the  maximum  power  density  is  103mWcm-2. 
The  four  curves  presented  in  Fig.  3  put  in  evidence  that  the  lim¬ 
iting  phenomenon  is  dry-out,  which  will  be  discussed  in  details 
in  the  following  section.  Both  the  model  and  experimental  results 
show  that  the  decrease  of  relative  humidity  at  high  current  den¬ 
sity,  which  causes  ionic  resistance  to  increase,  coincides  with  the 
non-linear  voltage  drop  that  stops  the  reaction. 

The  main  difference  between  experimental  and  analytical 
results  occurs  in  temperature  results.  However,  trends  linking 
temperature,  relative  humidity,  membrane  ionic  resistance  and 
polarization  curve  are  well  captured.  Considering  the  scope  of  the 
model,  which  was  to  demonstrate  the  manifestation  of  the  prob¬ 
lems  limiting  small  PEMFC  performance,  the  model  is  adequate. 

As  presented  in  Table  1,  fitting  parameters  were  necessary  to 
match  the  model  with  the  experimental  results.  /0  is  used  to  fit 
the  voltage  drop  from  reaction  kinetics  and  Cheat  is  a  heat  trans¬ 
fer  correction  factor.  Temperature  is  a  commonly  adjusted  value  in 
simplify  model.  The  heat  transfer  correlations  used  here  are  for  a 
highly  simplified  geometry,  neglect  radiation  and  density  changes 
from  concentration  gradients.  Thus,  they  are  certainly  not  quanti¬ 
tatively  exact  and  need  to  be  corrected  using  the  Cheat  factor. 

Based  on  the  modeling  results,  it  is  possible  to  estimate  the  rel¬ 
ative  importance  of  heat  loss  through  the  system.  The  mean  value 
of  the  ratio  between  Qc  and  Qc  +  Qout  is  10.7%.  This  small  fraction  of 


heat  loss  through  the  GDL  puts  in  evidence  the  importance  of  the 
packaging  on  heat  transfer.  The  fact  that  dry  hydrogen  was  supplied 
to  the  anode  using  a  stoichiometry  ratio  of  3  has  a  small  impact  on 
the  thermal  balance  since  heat  loss  through  heating  of  the  hydro¬ 
gen  supply  at  the  anode  is  only  up  to  1.1%  of  total  heat  loss  of  the 
system  and  heat  generated  through  condensation  of  water  is  less 
than  0.25%. 

4.  Investigation  of  water  management  problems 

4.1.  Dry-out 

As  presented  in  [13,17],  operation  of  a  small  PEMFC  is  divided 
into  three  main  regions:  humidification,  flooding,  and  dry-out. 
Fig.  4  presents  the  relative  humidity  curve  of  the  cathode,  which 
increases  up  to  100%  RHC  at  a  current  density  of  65  mA  cm-2 
(humidification).  In  the  second  region,  relative  humidity  is  sat¬ 
urated  at  100%,  which  means  that  water  is  condensing  in  the 
cell.  Water  flux  generated  through  the  electrochemical  reaction  is 
higher  than  water  flux  removed  via  free  convection.  At  a  current 
density  of  320  mA  cm-2  relative  humidity  of  the  cathode  begins 
to  decrease  with  increasing  current  density.  Under  these  condi¬ 
tions  the  increase  in  water  partial  pressure  with  current  density  is 
not  significant  enough  to  compensate  for  the  increase  in  saturation 
pressure  due  to  temperature  augmentation  of  the  cell. 

Moreover,  decreasing  relative  humidity  will  increase  ionic  resis¬ 
tance  of  the  electrolyte  and  lead  to  a  larger  voltage  drop  in  the  cell. 
This  loss  will  further  raise  the  fuel  cell  temperature  and  accelerate 
the  drop  in  relative  humidity  until  the  cell  reaches  a  limiting  cur¬ 
rent  (I\)  caused  by  an  excessively  high  ionic  resistance  and  the  heat 
generation  that  will.  This  third  zone  is  the  dry-out  of  the  cell  and 
has  been  explained  by  O’Hayre  et  al.  [17]. 

To  better  understand  the  manifestation  of  flooding  and  dry¬ 
out  phenomena,  Fig.  4  compares  the  water  flux  generated  by  the 
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Fig.  4.  Reference  model  results  for  water  flux  generated  by  electrochemical  reaction 
and  maximum  water  flux  removed  by  free  convection  mass  transfer.  Relative  humid¬ 
ity  at  the  cathode  is  also  plotted  as  a  function  of  current  density.  With  increasing 
current  density,  water  balance  in  the  cell  changes  from  humidification,  to  flooding, 
and  then  to  dry-out. 


electrochemical  reaction  versus  the  maximum  water  flux  removed 
through  free  convection  [17].  The  three  zones  in  Fig.  4  are  sepa¬ 
rated  by  two  points  of  operation  under  which  maximum  water 
removal  rate  is  at  the  same  level  as  the  water  generation  rate.  It 
should  be  noted  that  water  removal  flux  is  calculated  using  the 
saturated  value  of  water  in  air.  Therefore,  this  represents  the  maxi¬ 
mum  water  flux  that  can  be  removed  from  the  cell  by  evaporation. 
In  humidification  and  dry-out  regions,  where  water  flux  removal 
rate  at  100%  RH  is  higher  than  the  water  flux  generated,  the  water 
flux  removed  from  the  system  will  decrease  in  order  to  match  the 
water  flux  generated.  Under  this  condition,  relative  humidity  of  the 
cell  will  drop  under  100%  as  presented  in  relative  humidity  curve 


in  Fig.  4.  But,  in  the  flooding  region,  the  water  flux  removed  will 
not  be  able  to  match  the  generated  water  flux  because  it  is  already 
at  its  maximum  value.  Consequently,  steady-state  operation  in  this 
region  is  not  possible  without  flooding. 

The  transition  points  between  these  zones  are  dependant  on 
the  system  configuration  and  environmental  conditions.  Electri¬ 
cal  resistance,  thermal  resistance,  mass  transport  resistance,  MEA 
properties,  GDL  properties,  etc.  will  influence  the  position  and  man¬ 
ifestation  of  these  zones  of  operation.  In  addition  to  GDL  and  MEA 
properties,  packaging  of  the  system  will  have  a  strong  influence  on 
the  PEMFC  performance. 


4.2.  Thermal  resistance  variation 


Let  us  define  the  equivalent  thermal  resistance  (ftp)  of  a  small 
PEMFC  as  the  ratio  of  temperature  difference  between  the  cell  and 
the  environment  (AT)  over  heat  loss  from  the  cell,  averaged  over 
all  operation  points: 


Rj  = 


AT 

Qc  +  Qout 


(33) 


A  system  with  a  higher  thermal  resistance  will  have,  for  the 
same  amount  of  heat  generated  by  the  cell,  a  higher  temperature. 
With  results  presented  previously,  this  would  mean  that  increas¬ 
ing  the  thermal  resistance  between  the  PEMFC  and  its  environment 
would  increase  the  cell’s  temperature,  lower  relative  humidity  and 
ultimately  decrease  the  limiting  current  density  (I{)  by  pushing  for¬ 
ward  the  dry-out  problem.  This  point  has  already  been  highlighted 
by  O’EIayre  et  al.  [17]  through  modeling  analysis  while  consider¬ 
ing  the  effect  of  GDL  thickness.  But  in  a  more  comprehensive  way, 
this  phenomenon  not  only  occurs  while  changing  GDL  properties 
but  for  the  entire  packaging  of  the  system.  Lienee,  using  the  small 
PEMFC  prototype  presented  previously  as  a  reference,  tests  were 
run  while  changing  the  thermal  resistance  of  the  system  to  demon¬ 
strate  this  behaviour.  Experimental  results  for  the  same  cell  with 
the  surface  area  and  material  of  the  packaging  modified  to  lead  to 
higher  and  lower  thermal  resistance  are  presented  in  Fig.  5.  Ther- 


Fig.  5.  Experimental  results  for  free  convection  fuel  cell  with  different  thermal  resistance  values  (RT)  as  a  function  of  current  density,  showing  that  lowering  thermal  resistance 
increases  limiting  current:  (a)  polarization  curves,  (b)  relative  humidity  on  surface  of  cathode  GDL,  (c)  fuel  cell  membrane  resistance  calculated  by  Eqs.  (23)  and  (24),  based 
on  relative  humidity  at  the  surface  of  the  cathode  GDL,  and  (d)  temperature  on  surface  of  cathode  GDL  (ambient  air  at  23.5  °C  and  40%  RH). 
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mal  resistance  is  increased  by  using  material  with  lower  thermal 
conductivity  and  decreased  by  adding  conductive  material  to  add 
surface  area  in  contact  with  the  environment. 

As  shown  in  Fig.  5  decreasing  equivalent  thermal  resistance 
increases  the  limiting  current  (I\)  and  the  range  in  which  flood¬ 
ing  occurs.  The  increase  in  performance  is  only  due  to  delaying  the 
transition  point  between  the  flooding  and  dry-out  regions.  With 
higher  thermal  resistance,  temperature  variation  as  function  of  cur¬ 
rent  density  is  higher  and  the  saturation  pressure  increase  is  more 
significant.  Thus,  the  humidity  level  decreases  and  the  ionic  resis¬ 
tance  increases  are  more  significant  which  causes  dry-out  to  occur 
earlier.  This  behaviour  is  pictured  by  the  four  curves  presented  in 
Fig.  5. 

Using  the  model  developed,  simulations  of  cells  with  various 
equivalent  thermal  resistances  were  carried  out  to  obtain  the  lim¬ 
iting  current  density  (I\).  Changing  thermal  resistance  in  the  model 
slightly  affects  the  mass  transport  resistance  by  changing  the  tem¬ 
perature  gradient  which  affects  air  properties.  Fig.  6  shows  that 
the  increase  in  limiting  current  density  as  function  of  equivalent 
thermal  resistance  variation  is  highly  non-linear  and  increases  dra¬ 
matically  around  frr  =  25KW_1.  As  expected  from  experimental 
results,  lowering  thermal  resistance  increases  the  limiting  current 
up  to  a  maximum  limiting  current  density ,  /i_max  (675  mA  cm-2  in 
Fig.  6).  The  maximum  limiting  current  density  (/i_max)  is  the  current 
value  for  which  voltage  is  zero  and  the  limiting  water  manage¬ 
ment  phenomenon  is  at  the  point  of  transitioning  between  flooding 
and  dry-out.  At  this  operating  point,  a  small  decrease  in  the  equiv¬ 
alent  thermal  resistance  would  not  allow  the  system  to  heat  up 
sufficiently  to  reach  dry-out  operation  mode  before  voltage  drops 
to  OV.  Consequently,  under  these  conditions,  the  limiting  factor 
for  the  PEMFC  would  not  be  dry-out  but  flooding  instead.  Limit¬ 
ing  current  would  then  be  shifted  to  the  transition  point  between 
the  humidification  and  flooding  zones  in  Fig.  4  (crossing  point 
furthest  to  the  left).  The  important  point  shown  in  Fig.  6  is  that 
thermal  resistance  can  significantly  increase  the  limiting  current  of 
the  cell  and  control  whether  the  limiting  phenomenon  is  dry-out  or 
flooding. 

4.2.1.  Low  thermal  resistance  behaviour 

If  thermal  resistance  of  the  system  is  low  enough,  increasing 
heat  loss  will  not  cause  the  temperature  of  the  cell  to  rise  sig¬ 
nificantly  and  consequently  the  saturation  pressure  will  remain 
practically  constant.  Increasing  the  current  density  will,  however, 


Fig.  6.  Model  results  showing  the  impact  of  equivalent  thermal  resistance  (RT)  on 
the  limiting  current  density  for  free  convection  fuel  cells.  As  the  RT  is  increased,  the 
limiting  phenomenon  switches  from  flooding  to  dry-out  (ambient  air  at  20-30 °C 
and  40%  RH). 


increase  rate  of  water  generated  and  the  water  partial  pressure.  At 
practically  constant  temperature,  water  will  eventually  accumulate 
in  the  system  without  manifestation  of  the  dry-out  phenomenon. 
Experimentally,  this  operating  mode  can  be  achieved  by  putting 
the  anode  side  of  the  small  PEMFC  in  contact  with  a  large  heat  sink, 
such  as  a  massive  block  of  thermally  conductive  material  (stainless 
steel).  This  operating  mode  may  seem  unusual  for  a  small  PEMFC, 
but  in  fact  it  simulates  a  cell  that  would  be  physically  integrated 
into  a  large  system.  This  effect  would  be  even  greater  for  micro  fuel 
cells. 

The  experimental  results  presented  in  Fig.  7  show  that  tem¬ 
perature  variation  on  the  surface  of  the  GDL  of  a  small  PEMFC  in 
contact  with  a  heat  sink  is  very  low;  temperature  variation  is  less 
than  0.5  °C  for  a  current  density  up  to  100  mA cm-2.  Under  these 
operating  conditions,  the  limiting  current  in  steady-state  operation 
is  low  (150  mA  cm-2 )  and  occurs  because  of  the  formation  of  water 
droplets  on  the  surface  of  the  GDL  as  observed  visually;  the  dry¬ 
out  phenomenon  is  never  reached.  This  highlights  the  importance 
of  thermal  rejection  on  water  management  problems,  which  in  this 
case  creates  a  flooding  problem  instead  of  a  dry-out  phenomenon. 

4.2.2.  Robustness 

As  shown  in  Fig.  6,  change  in  environment  temperature  or  rel¬ 
ative  humidity  shifts  the  limiting  current  curve  and  changes  the 
optimal  point  of  equivalent  thermal  resistance.  In  order  to  achieve 
system  robustness,  one  would  have  to  consider  these  phenomena 
and  use  a  more  indulgent  heat  rejection  configuration  in  order 
to  achieve  dry-out  phenomenon  on  a  larger  range  of  values  and 
avoid  flooding  as  the  limiting  phenomenon.  Otherwise  a  system 
designed  tightly  around  an  optimum  value  for  a  certain  environ¬ 
mental  condition  might  not  be  able  to  reach  steady-state  operation 
in  environments  of  higher  relative  humidity  and/or  lower  temper¬ 
ature. 


4.3.  Mass  transport  resistance  variation 


Let  us  define  the  equivalent  mass  transport  resistance  ( Rm  )  of  a 
small  PEMFC  as  the  ratio  between  water  partial  density  variation 
and  water  generated  by  the  system,  averaged  over  all  operation 
points: 


Rm  = 


o 


mH7o 


(34) 


As  for  thermal  resistance,  mass  transport  resistance  will  impact 
the  performance  of  a  small  PEMFC  running  in  free  convection  mode. 
If  system  mass  transport  resistance  is  increased,  the  water  molar 
fraction  will  increase  and  the  oxygen  molar  fraction  at  the  cath¬ 
ode  will  decrease.  From  a  water  management  perspective,  this  will 
allow  the  cell  to  reach  higher  current  density  due  to  lower  water 
removal  rate,  which  will  push  forward  the  transition  point  between 
flooding  and  dry-out.  This  will  be  true  unless  mass  transport  resis¬ 
tance  is  too  high  and  creates  significant  voltage  drop  due  to  lack  of 
oxygen.  Using  the  modeling  results,  we  can  plot  oxygen  and  water 
vapour  molar  fraction  against  current  density. 

Fig.  8  shows  that  the  oxygen  molar  fraction  does  not  decrease 
significantly,  even  at  high  current  density.  Molar  fraction  of  oxy¬ 
gen  varies  from  0.21  to  0.14  when  current  density  goes  from  0 
to  300  mA cm-2.  Also,  according  to  this  trend,  the  model  predicts 
that  oxygen  molar  fraction  would  reach  0  at  about  1600  mAcnrr2, 
which  is  quite  superior  to  any  results  presented  here  and  for  any 
free  convection  micro-PEMFC.  An  important  conclusion  that  can  be 
drawn  is  that  oxygen  depletion  at  the  cathode  catalyst  zone,  due  to 
pure  mass  transport  problems,  is  not  the  actual  problem  that  lim¬ 
its  small  PEMFC  when  running  in  free  convection  mode.  The  same 
conclusion  was  drawn  by  Modroukas  [15]  and  O’Hayre  et  al.  [17]. 
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Current  density  [mA/cm2] 


Fig.  7.  Experimental  results  for  free  convection  fuel  cell  operated  with  a  high  thermal  resistance  and  a  low  thermal  resistance  (heat  sink)  that  allows  practically  constant 
temperature  as  a  function  of  current  density,  showing  that  low  thermal  resistance  system  was  limited  by  flooding  problems:  (a)  polarization  curves,  (b)  relative  humidity  on 
surface  of  cathode  GDL,  (c)  fuel  cell  membrane  resistance  calculated  by  Eqs.  (23)  and  (24),  based  on  relative  humidity  at  the  surface  of  the  cathode  GDL,  and  (d)  temperature 
on  surface  of  cathode  GDL  (ambient  air  at  23.5  °C  and  40%  RH). 


The  experimental  impact  of  changing  mass  transport  resistance 
was  evaluated  by  using  as  reference  the  small  PEMFC  presented  pre¬ 
viously  and  by  adding  a  screen  on  the  cathode’s  GDL  top  surface. 
The  screen,  which  consists  of  a  carbon  cloth  identical  to  the  cur¬ 
rent  GDL,  increases  the  mass  transport  resistance  to  and  from  the 
cathode.  For  these  tests,  temperature  and  relative  humidity  were 
measured  at  the  top  surface  of  the  cathode  GDL  in  order  to  compare 
with  previous  results. 


Fig.  8.  Reference  model  results  for  relative  humidity  and  species  molar  fraction  at 
the  surface  of  cathode  GDL  as  a  function  of  current  density  for  free  convection  fuel 
cell  operated  at  ambient  air  at  23.5  °C  and  40%  RH.  These  results  show  that  oxygen 
depletion  is  not  the  limiting  phenomenon,  even  at  higher  current  densities. 


The  trends  described  above  are  shown  experimentally  in  Fig.  9. 
Adding  mass  transport  resistance  to  the  system  increases  maxi¬ 
mum  limiting  current.  With  a  higher  mass  transport  resistance  the 
system  will  have  a  higher  molar  fraction  of  water  at  the  cathode, 
so  for  a  same  temperature,  the  water  activity  in  the  membrane  will 
be  higher.  This  will  prevent  dry-out  of  the  electrolyte  and  increase 
the  limiting  current.  Fig.  9  shows  that  at  285  mA cm-2,  the  rela¬ 
tive  humidity  of  the  reference  test  is  decreasing  while  the  relative 
humidity  of  the  prototype  with  higher  mass  transport  resistance 
is  still  saturated  at  100%.  Up  to  this  point  the  temperature  val¬ 
ues  for  both  tests  are  similar  but  because  of  the  relative  humidity 
decrease,  the  reference  test  temperature  starts  to  increase  due  to 
augmentation  of  ionic  resistance  that  generates  more  heat  loss. 
Adding  a  screen  on  the  surface  of  the  GDL  increases  the  thermal 
resistance,  and  as  presented  above,  this  should  decrease  the  limit¬ 
ing  current.  But,  the  increase  in  mass  transport  resistance  is  more 
significant  and  performances  increase  instead  of  decreasing.  This 
proves  that  increasing  mass  transport  resistance  of  a  small  PEMFC 
that  is  nominally  limited  by  dry-out  can  increase  its  performances 
by  keeping  water  in  the  cell.  Just  as  thermal  resistance  was  found 
to  be  critical,  changing  mass  transport  resistance  can  also  mod¬ 
ify  the  fuel  cell  performance  and  impact  the  water  management 
challenges. 


5.  Discussion  on  water  management 

5.1.  Limiting  phenomena 

Reducing  thermal  resistance  and  increasing  mass  transport 
resistance  allows  small  air  breathing  PEMFC  to  reach  higher  lim¬ 
iting  current  densities.  This  is  done  by  pushing  the  transition 
point  between  flooding  and  dry-out  phenomena  to  higher  current 
densities.  However,  the  operating  range  in  which  flooding  occurs 
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Fig.  9.  Experimental  results  for  free  convection  fuel  cell  with  different  mass  transport  resistance  (RM),  showing  that  increasing  mass  transport  resistance  increases  the 
limiting  current:  (a)  polarization  curves,  (b)  relative  humidity  on  surface  of  cathode  GDL  as  a  function  of  current  density,  (c)  fuel  cell  membrane  resistance  calculated  by  Eqs. 
(23)  and  (24),  based  on  relative  humidity  at  the  surface  of  the  cathode  GDL  as  a  function  of  current  density,  and  (d)  temperature  on  surface  of  cathode  GDL  as  a  function  of 
current  density  (ambient  air  at  23.5  °C  and  40%  RH). 


broadens.  Ultimately,  with  further  reduction  in  thermal  resistance, 
dry-out  operating  mode  will  not  be  reachable  and  only  flooding 
of  the  system  will  be  achievable  beyond  the  humidification  oper¬ 
ating  zone.  The  packaging,  through  its  thermal  and  mass  transport 
resistance,  establishes  what  kind  of  limiting  phenomena  will  occur. 
For  example,  a  small  package  with  high  thermal  resistance  and  low 
mass  transport  resistance  will  exhibit  the  dry-out  phenomenon. 
In  the  opposite  situation,  a  large  package  with  low  thermal  resis¬ 
tance  and  high  mass  transport  resistance  can  result  in  flooding 
problems.  This  conclusion  allows  to  unify  observations  from  Refs. 
[14-16]  stating  that  the  limiting  phenomenon  is  flooding  with 
conclusions  from  Refs.  [8,13,17]  that  the  limiting  phenomenon  is 
dry-out. 

However,  another  point  that  must  be  considered  regarding  the 
manifestation  of  limiting  phenomena  is  the  ideal  limiting  current 
density  (/i_ideai)»  which  is  the  limiting  current  of  the  system  if  no 
mass  transport,  flooding  or  dry-out  phenomena  occur.  The  limiting 
current  value  in  this  ideal  case  is  solely  due  to  the  voltage  drop  from 
reaction  kinetics  and  Ohm’s  law  with  constant  ionic  resistance,  so 
this  value  is  fixed  by  cell  properties  and  electrical  resistance  of  the 
system.  This  value  will  become  the  limiting  phenomenon  when  the 
following  ratio  is  respected: 


L^eal  <  1  (35) 


In  this  particular  case,  the  limiting  current  (I\)  due  to  dry-out  or 
flooding  will  not  be  reachable  because  loss  from  electrical  resis¬ 
tances  or  reaction  kinetics  would  be  too  important.  In  order  to  get 
higher  performance,  one  would  have  too  lower  these  losses  instead 
of  working  on  thermal  or  mass  transport  resistance. 


5.2.  Operating  considerations 

Using  the  model,  relative  humidity  for  three  different  thermal 
resistances  was  plotted  in  Fig.  10.  Sufficiently  low  thermal  resis¬ 
tance  will  maintain  the  flooding  operating  mode  without  reaching 
the  dry-out  mode,  as  shown  for  RT  =  11.6  in  Fig.  10.  This  condition 


Fig.  10.  Model  results  of  relative  humidity  on  surface  of  cathode  GDL  as  a  func¬ 
tion  of  current  density  for  free  convection  fuel  cell.  For  high  equivalent  thermal 
resistance  (RT  =  59.1 1<  W_1 )  there  is  no  flooding  zone,  only  humidification  and  dry¬ 
out  occur.  For  intermediate  equivalent  thermal  resistance  (RT  =  41.5KW_1)  there 
are  three  zones:  humidification,  flooding  and  dry-out.  For  low  equivalent  thermal 
resistance  (RT  =  11.6  K  W_1 )  there  is  no  dry-out  zone  (ambient  air  at  23.5  °C  and  40% 
RH). 
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is  possible  if  temperature  variations  are  kept  very  low  so  that  the 
partial  pressure  of  water  will  remain  higher  than  the  saturation 
pressure  value.  Under  this  condition,  the  transition  point  between 
flooding  and  dry-out  operating  mode  disappears  and  the  limiting 
current  of  the  system  is  fixed  by  the  transition  between  humidifi¬ 
cation  and  flooding  operating  mode.  The  system  could  then  only  be 
operated  in  steady-state  regime  in  the  humidification  zone. 

However,  if  thermal  resistance  is  increased  sufficiently,  the 
flooding  operating  zone  will  disappear,  as  presented  for  the  equiva¬ 
lent  thermal  resistance  curve  of  Rj  =  59.1  K  W-1  in  Fig.  10.  Maximum 
water  removal  rate  will  always  be  greater  than  water  generated  by 
the  electrochemical  reaction.  Water  flux  removed  will  adjust  itself 
to  the  water  generated  and  will  bring  the  system  down  to  a  relative 
humidity  lower  than  100%.  Considering  the  shape  of  the  maximum 
water  removal  rate  at  low  current  density,  relative  humidity  will 
first  increase  without  reaching  100%  and  eventually  will  decrease. 
Therefore,  only  humidification  and  dry-out  operating  modes  will 
occur.  However,  it  is  important  to  consider  that  every  operating 
point  will  be  achievable  because  there  is  no  flooding  zone,  which 
is  not  the  case  with  the  intermediate  equivalent  thermal  resistance 
curve  (Rj  =  41.5  KW-1).  Under  this  mode  of  operation,  there  will 
be  a  range  in  which  long-term  operation  is  not  possible  because 
of  flooding.  In  order  to  decrease  the  flooding  zone  and  increase 
robustness,  the  thermal  resistance  must  be  increased  and/or  the 
mass  transport  resistance  must  be  decreased.  However,  increasing 
robustness  of  the  system  is  done  at  the  cost  of  reducing  the  limiting 
current  density. 

5.3.  Thermal  resistance  and  mass  transport  resistance  ratio 

The  previous  discussion  for  water  management  within  free 
convection  fuel  cells  was  mainly  conducted  considering  various 
thermal  resistances  but  at  a  constant  mass  transport  resistance. 
However,  as  demonstrated  similar  results  can  be  achieved  by 
changing  mass  transport  resistance  and  keeping  the  thermal  resis¬ 
tance  practically  constant.  Hence,  increasing  thermal  resistance 
and/or  lowering  mass  transport  resistance  reduces  the  humidity 
level  in  the  system  either  by  increasing  the  saturation  pressure  or 
lowering  the  water  partial  pressure.  This  favours  dry-out  as  the 
limiting  phenomenon.  On  the  opposite,  decreasing  thermal  resis¬ 
tance  and/or  increasing  mass  transport  resistance  would  lead  to  a 
flooding  problem.  Since  these  resistances  variations  act  in  opposite 
ways,  one  can  define  a  ratio  between  thermal  and  mass  transport 
resistances  that  would  describe  which  limiting  phenomenon  will 
occur.  Considering  the  units  of  thermal  and  mass  transport  resis¬ 
tance,  we  can  define  a  non-dimensional  water  management  ratio  of 
the  system  (WMratio)  by  multiplying  the  ratio  of  resistances  by  the 
specific  heat  and  density  of  air: 

WMrati0  =  ^CPP  (36) 

For  the  reference  prototype  presented  in  this  paper,  the  transition 
point  between  dry-out  and  flooding  limitation  phenomena  would 
occur  at  WMratio  =  0.026  for  a  PEMFC  operated  at  ambient  air  of 
23.5  °C  and  a  relative  humidity  of  40%.  It  is  important  to  understand 
that  this  specific  value  is  valid  only  for  the  prototype  and  its  specific 
geometry. 

6.  Conclusion 

In  this  paper,  we  have  presented  a  simple  ID  model  of  a 
free  convection  PEMFC  that  couples  heat  and  mass  transfer  with 
the  electrochemical  reaction  of  a  small  PEMFC.  Experimental 
results  using  a  small  air  breathing  PEMFC  prototype  with  a  planar 
configuration  validate  the  consistency  of  this  model.  However,  sim¬ 
plifications  used  in  the  model  regarding  heat  transfer  did  not  allow 


an  exact  agreement  between  experimental  and  model  results,  so 
a  correction  factor  was  necessary  to  get  a  better  match  between 
temperatures.  The  model  would  have  advantage  to  be  improved 
regarding  heat  transfer  analysis.  Also,  the  analogy  between  heat 
and  mass  transfer  was  used  to  calculate  coefficients  and  each 
phenomenon  was  considered  separately.  The  model  would  be 
improved  if  one  would  merge  these  two  phenomena  while  eval¬ 
uating  heat  and  mass  transfer  coefficients. 

Analysis  was  conducted  using  the  model  along  with  various 
experimental  investigations  to  evaluate  dry-out  and  flooding  phe¬ 
nomena  and  the  impact  of  thermal  resistance  and  mass  transport 
resistance  of  a  small  PEMFC  running  in  free  convection  mode. 
Discussion  of  water  management  was  made  regarding  limiting  phe¬ 
nomena  and  operation  of  such  a  system.  Dependency  of  thermal 
resistance  and  mass  transport  resistance  of  the  system  was  demon¬ 
strated  concerning  cell  performances  limitations  due  to  water 
management  problems. 

Important  conclusions  were  drawn  concerning  water  manage¬ 
ment  problems  for  small  free  convection  fuel  cells: 

(1) The  ratio  between  thermal  and  mass  transport  resistance 
(WMratio)  of  the  packaging  is  a  critical  aspects  for  water  man¬ 
agement,  controlling  if  and  when  dry-out  or  flooding  problems 
will  occur. 

(2 )  High  ratio  of  resistances  tends  to  dry-out  the  electrolyte,  leading 
to  increased  ionic  resistance  and  a  voltage  drop  limiting  the 
achievable  current  density. 

(3)  Low  ratio  of  resistances  could  lead  to  flooding  of  the  cathode 
that  would  prevent  oxygen  mass  transport  to  the  catalytic  sites 
and  limit  the  achievable  current  density.  This  is  characteristic 
of  small  and  micro  fuel  cells  integrated  into  larger  system. 

(4)  Fuel  cells  running  in  free  convection  mode  are  characterized 
by  three  modes  of  operation:  humidification,  flooding  and  dry¬ 
out,  as  presented  by  Refs.  [13,17].  The  manifestation  of  these 
modes  is  affected  by  the  system’s  packaging  properties  and  the 
surrounding  conditions. 

(5)  Considering  fuel  cells  limited  by  dry-out,  there  is  a  compromise 
between  the  maximum  limiting  current  density  and  robustness 
of  the  system:  the  flooding  zone  will  broaden  when  the  maxi¬ 
mum  current  density  is  increased  through  a  diminution  of  the 
ratio  between  thermal  and  mass  transport  resistance. 

(6)  Mass  transport  voltage  drop  caused  by  a  lack  of  oxygen  at  the 
cathode,  due  to  pure  mass  transport  problems  and  not  flooding, 
is  not  significant  considering  the  limiting  current  density  of  air 
breathing  small  PEMFC.  The  same  conclusion  was  also  drawn 
by  Modroukas  [15]  and  O’Hayre  et  al.  [17]. 

Finally,  the  conclusions  brought  by  this  paper  show  the  impact 
of  packaging  for  performances  reached  by  a  small  fuel  cell.  The 
importance  of  heat  and  mass  transfer  of  an  air  breathing  fuel  cell 
makes  the  design  of  the  packaging  critical  for  the  performance  of 
the  system.  Thus,  further  work  concerning  small  PEMFC  would  ben¬ 
efit  from  placing  emphasis  on  packaging  design.  This  is  even  more 
critical  for  micro  fuel  cells  which  are  dominated  by  heat  transfer 
through  the  packaging.  For  low  thermal  resistance  systems  empha¬ 
sis  should  be  placed  on  removing  adequately  water  in  order  to 
increase  performance.  On  the  other  hand,  for  high  thermal  resis¬ 
tance  system,  emphasis  should  be  placed  on  increasing  thermal 
rejection  and  lowering  water  removal  to  achieve  higher  perfor¬ 
mance  [17].  One  would  need  to  keep  in  mind  that  thermal  rejection 
and  mass  transport  are  linked  by  the  free  convection  phenomena 
which  will  push  both  elements  in  the  same  direction.  However,  it 
is  possible  to  change  more  significantly  one  element.  For  example 
mass  transport  resistance  will  be  more  affected  than  thermal  resis¬ 
tance  by  changes  in  GDL  porosity,  while  changing  the  packaging  size 
without  changing  the  cell’s  active  surface  or  cathode  exposed  area 
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will  affect  the  thermal  resistance  more  considerably.  The  simple 
modeling  approach  and  non-dimensional  water  management  ratio 
defined  here  can  guide  a  designer  in  developing  high  performance 
air  breathing  fuel  cells. 
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